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Abstract. ISO spectral observations of the supernova 
remnant RCW103 are presented. This object is the pro- 
totype of relatively young remnants (~ 10 3 yr) with fast 
shocks (v s ~ 1000 km/s) interacting with dense interstellar 
medium. 

The spectrum is dominated by prominent lines of 
[Nell] , [Sill] , [Fell] and other low excitation species which 
provide, for the first time, a simple and reliable estimate 
of the gas abundances of refractory (Si, Fe, Ni) and non- 
refractory (Ne, P, S, CI, Ar) species. Apart from nickel, 
all the derived abundances are close to solar, confirming 
that the shock has destroyed all dust grains. Like the op- 
tical nickel lines, [NiII]A6.64 /im yields Ni abundances a 
factor ~10 solar which we propose results from a large 
underestimation of the computed Ni + collision strengths. 

The observed intensities and velocity widths of ionic 
lines are compatible with emission from the post-shock 
region alone with only a very small (if any) contribu- 
tion from the photoionizcd precursor. This result does not 
agree with shock models which predict that the precursor 
should emit powerful line emission, especially from highly 
ionized species. The possible consequence of this on the 
modelling of Seyfert spectra is briefly discussed. 



Key words: Atomic data; ISM: abundances; ISM: su- 
pernova remnants; Galaxies: Seyfert; Infrared: ISM: lines 
and bands 



1. Introduction 

Radiative supernova remnants are the ideal astrophysi- 
cal laboratory for studying the emission spectrum of rela- 
tively fast shocks interacting with the interstellar medium 



(e.g. Draine & McKee 1993). In the standard, idealized, 
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view the supernova blast wave expands through a uniform, 
low density (~1 cm -3 ) medium and the shock remains 
adiabatic, i.e. the postshock gas has insufficient time to 
cool/recombine, for > 10 4 yr and up until the shock has 
slowed to v s < 200 km/s. After this the shock becomes 
radiative with the postshock region recombining and ra- 
diating away most of the shock mechanical energy in the 
form of UV, optical and IR lines whose surface bright- 
ness simply scale with the shock energy flux nv 3 , n be- 
ing the pre-shock density. In practice, a 'standard' radia- 
tive supernova remnant is expected to emit weak lines, i.e. 
E(H/3)<10^ 5 erg cm~ 2 s" 1 sr" 1 (cf. Sect 3.2 of Dopita & 
Sutherland 1996, hereafter DS96) and 2-3 orders of mag- 
nitude weaker than those observed in bright supernovae 
remnants such as RCW103 (Oliva et al. 1989, hereafter 
pMD89| ). This simple fact indicates, therefore, that the 
shock of RCW103 is much faster or, equivalently, that the 
SNR is much younger than the canonical values. This is 
also confirmed by X-ray spectra which yield an age of 
only 10 3 yr and a main shock velocity of about 1200 km/s 
(Nugent et al. |1984|). 



Thus, the prominent line emission from the optical fila- 
ments of RCW103 results from the impact of the v s ~ 1000 
km/s SNR blast wave onto dense clouds of ISM material 
(probably molecular clouds). Due to the higher density, 
the postshock gas cools and recombines much faster than 
in the canonical model and the SNR becomes radiative 
at much earlier times. The large surface brightness of the 
emitted line spectrum simply reflects the large mechanical 
energy flux of the shock. 

Studying the IR spectrum of radiative shocks can give 
interesting information on the metal abundances and, by 
comparing refractory and non-refractory species, can yield 
a direct measurement of the actual efficiency of grain de- 
struction by the shock. Moreover, the spectra can be used 
to verify the predictions of mod els of fast shocks, which 
have been recently published by DS96. 
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Fig. 1. Individual line scans with SWS02. Wavelengths are in and fluxes in Jy. Arrows mark the positions of 
undetected lines. 



Being the brightest and best studied SNR, RCW103 
was chosen as the prototype object for ISO template spec- 
tra of gas excited by fast shocks. The observations are de- 
scribed in Sect. ^ and the derived physical parameters are 
discussed in Sect. [j[ The results are compared with shock 
model predictions in Sect. ^ where we draw particular at- 
tention to the relative roles of the photoionized precursor 
and post-shock region in producing the observed lines. 



2. Observations and results 

2.1. ISO spectroscopy 

Complete SWS01 (speed 4, 6700 sec total integration 
time) and LWS01 (4400 sec total integration time, no 
background spectrum subtracted) spectra centered at the 
peak of [FcII]1.644 /j,m line emission (cf. Fig. ^) were 
obtained on February 6 and February 20, 1996. These 
were complemented by a quick PHT-S spectrum (448 sec, 
February 20) and deeper SWS02 observations at selected 
wavelengths (total 7000 sec, obtained on August 15, 1996 
and February 17, 1997), always centered at the same po- 
sition. The short wavelength section (2.5-5 /im) of the 
PHT-S spectrum was strongly contaminated by detector 
memory effects (a very bright source was observed just 
before RCW103), the PHT-SS results are therefore unre- 
liable and not presented here. 

The SWS data were reduced using standard routines 
of the SWS interactive analysis system (IA) using calibra- 
tion tables as of September 1997. Reduction relied mainly 
on the default pipeline steps, plus removal of signal spikes, 
elimination of the most noisy band 3 detectors, and flat- 
fielding. The LWS spectrum is based on the end-product 



of the automatic pipeline as of April 1997 (i.e. OLP 6). 
A post-processing was performed within the ISO Spec- 
tral Analysis Package (ISAP)|^|, Version 1.2, with special 
emphasis on removal of signal spikes and memory effects, 
averaging of the different scans, and flat-fielding of the 10 
detectors. 

The final rebinned spectra are displayed in Figs, [l] to 
[| and the derived line fluxes are listed in Table [| together 
with additional information. Note that most of the lines in 
the 2.4-40 /im range were observed twice, i.e. in the com- 
plete SWS01 spectrum (Fig. |) and in the deeper SWS02 
line scans (Fig. [l]) . The derived fluxes were always within 
±30% and in most cases equal to much better than 20%. 
The errors quoted in Table [l] also include the differences 
between the two measurements. It should be noted that, to 
the best of our knowledge, the transitions of [C1II]A14.36 
and [PII]A32.8 are newly detected astronomical lines. Of 
interest is also the marginal detection of an unidentified 
feature at 74.2 /im whose position and flux are remarkably 
similar to those found in the spectrum of NGC7027 (Liu 
et al. |1996| ). 

The SWS spectrum (Figs, [|, ||) is characterized by 
prominent lines over a faint continuum (about 0.5 Jy at 
10 /itm) while emission by cold dust is evident in the LWS 
spectrum (Fig. ||). The level of the 100 /im continuum is 
similar to the background IRAS level reported by Arendt 
( 1989| ) and the continuum seen in the LWS spectrum is, 
therefore, probably dominated by back/foreground emis- 



The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 
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Fig. 5. Line images at various velocity bins, reconstructed from the IRSPEC spectra discussed in Sect. |[ all frames 
are cut to the same levels as the v=-60 km/s images. Note that [Fell] (1.644 /im) shows bright filaments moving by 
up to ±250 km/s, while H2 (2.121 fim) is narro w and unresolved (i.e. FWHM<130 km/s). The colour images are also 
available at http: / /www. arcetri.astro.it / ^oliva 



sion from the galaxy disk. The same applies to the PAH's 
features visible in the PHT-S spectrum (Fig. |J) and to the 
[OIJA63.2, [NII]A121 and [CII]A157 lines which are likely 
to be strongly contaminated by emission from the diffuse 
ISM. 



Given the relatively large extinction towards RCW103, 
i.e. A H ~ 0.75 or A v ^4.5 (cf. Oliva et al. 1990, hereafter 



OMD90 ), the reddening corrections are not negligible for 
the lines at the shortest wavelengths. Therefore, in Ta- 
ble |] we also list the extinctions values which were derived 
assuming a 'typical' reddening curve, i.e. rwA^ 1 ' 75 out- 
side of the silicate band. The largest and most uncertain 
correction is for [SIV]A10.52 which lies within the silicate 
band and for which we have assumed A(10/im)/Av— 0.1. 



Useful dynamical information can be derived from the 
profiles of the lines between 15 and 19 //m, the wavelength 
range where SWS02 grating spectra achieve their highest 
spectral resolution. The instrumental resolution depends 
on the size of the source along the dispersion direction and 
varies, for example, between 120 and 175 km/s for a com- 
pact and extended [SHI] A 18. 7 source, respectively. Luck- 
ily, the SWS slit was almost exactly aligned N-S and was, 
therefore, roughly uniformly illuminated in the dispersion 
(E-W) direction (cf. Fig. [?]). The observed line profiles 
arc displayed in Fig. [| The ionic lines are resolved and 
exhibit similar profiles, within the errors, but are broader 
than the H2 line which is unresolved. This agrees well with 
the higher spatial resolution velocity maps from the NIR 
data discussed below. 



2.2. Ground based near infrared spectroscopy 



Near infrared observations were collected in Jan 1992 at 
the ESO-NTT using the long-slit spectrometer IRSPEC 
equipped with a 62x58 SBRC InSb array whose pixel size 
was 2.2" along the slit and ~5 A along the dispersion 
direction. Line images of [Fell] Al. 644 and H2 A2.121 were 
reconstructed from 57 spectra with the 100" x 2.2" slit 
aligned E-W and shifted by steps of 2.2" along the N- 
S direction. Each long-slit spectrum consisted of a single 
on-chip integration of 30 sec with sky exposures every 10 
spectra. 

The integrated line images are displayed in Fig. [?] 
whose caption also include results of [Fell] aperture pho- 
tometry which has been used to determine the correction 
factors for the different beams used by SWS (see also Ta- 
ble |Tj) . These assume that all ionized species have spatial 
distribution similar to [Fell] which is justified by the fol- 
lowing arguments. 

- The flux of Bra seen by ISO is within <10% of the 
value extrapolated from the ground based measurement 
( PMD90Q assuming a constant [Fell] /Bra ratio over the 
region of interest. 

- The morphology of the [Fell] filaments is virtually iden- 
tical to those seen in optical line images (cf. e.g. Moorwood 
et al. |1987| ). 

It should be noted, however, that the H2 lines arise from 
a totally different region ~20" south of the ionized gas 
and outside the optical/radio/X-ray remnant, as origi- 



nally found by OMD89 This indicates that H2 emission 



traces material which has not yet been reached by the 
shock, most probably a molecular c loud hea ted by the 
soft X-rays from the shock front (cf. |OMD90| ) . 



Oliva et al.: ISO infrared spectroscopy of fast shocks. I. Ionized lines in RCW103 



Table 1. Observed ISO line fluxes 
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(2) Extinction (mag) extrapolated from A H =0.75 ( pMD90[ ) 
using Aa~A~ 1,7 up to 9 /xm. Reddening is assumed to be 
negligible beyond 12 [im 

(3) Size (in arcsec) of the SWS (A<45 /xm) and LWS (A>45 
/im) apertures 

Correction factor to account for different apertures, based 
on the IRSPEC map of [FeII]A1.644 (see caption of Fig. Q). 
a Only SWS01 measurement available 

Adopting a "standard" A(silicate feature) /Av — 0.1 ratio 

c Line flux is probably contamined by back/foreground 
galactic emission 

d Unidentified feature also seen in NGC7027 (Liu et al. |l996|) 
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Fig. 2. Complete 2.4-45 /im SWS01 spectrum with line 
identification. Note that each spectral segment has been 
"continuum subtracted" to remove instrumental offsets 
and drifts. The continuum level can be estimated from 
the PHT-S spectrum (Fig. |) and from the deeper SWS02 
line scans (Fig. |l|). 

Fig. shows images at various velocity bins, each 
roughly corresponding to the wavelength range covered 
by 1 pixel. Evident are the high velocity [Fell] filaments 
whose projected velocities extend up to ±250 km/s and 
are compatible with the idea that this line is produced 
downstream of the fast shock (cf. the Introduction). The 
H2 filaments, on the contrary, do not show evidence of 
motions larger than the FWHM=130 km/s instrumental 
resolution. 



3. Physical parameters of the emitting gas 

3.1. Temperature and density 

The ISO spectra include several density and/or temper- 
ature sensitive line pairs. Besides the well known, n e - 
sensitive [OIII], [Nelll], [NeV] and [SIII] fine structure 
lines, useful transitions are those of [Fell] whose paths are 
displayed in Fig. |q, while their dependence on T e and n e is 
visualized in Fig. Note that the ratios [FeII]A35.3/A26.0 
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tion. Note that RCW103 lies close to the the Galactic 
plane (b = —0.36). Therefore, the 100 fim continuum and 
the [OI], [CII], [Nil] lines are strongly contaminated, and 
most probably dominated by fore/background emission. 
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Fig. 4. Complete PHT-SL spectrum. The PAH features 
are most probably dominated by fore/background emis- 
sion from the Galaxy disk. 



and [FeII]A24.5/A17.9 are also density sensitive, but vary 
by only a factor ~1.7 between the low and the high density 
limits so the resulting n e is highly uncertain. 

The deduced values of T e and n e are summarized in 
Table || where the most remarkable result is the large 
spread of densities deduced from the various diagnostic 
ratios. The lowest value is from the [OIII] lines which 
are measured through the much larger LWS beam (cf. 
Table [l]) and are probably affected by the emission of a 
diffuse, lower density component. Some of the discrepan- 
cies between the [Fell] densities could also reflect density 
and temperature stratification. In particular, the ground 
state line at 26.0 /im could partially arise from excitation 
by collisions with atomic hydrogen in regions far down- 
stream of the shock where the temperature is too low 
to produce the other [Fell] lines. However, such a cold 
component should primarily show up in the temperature 
sensitive [Fell] A26.0/A1.64 ratio which, on the contrary, 
yields a relatively high value for T e and therefore implies 
that any cold component can account, at most, for 30% 




[Ne III]X15.56 
S(l) A17.04 
[Fe IIJA17.93 
[S III]X18.71 



-400 -200 200 400 

Velocity (km/s) 

Fig. 6. Normalized profiles of FIR lines falling at wave- 
lengths where the SWS spectrometer achieves its high- 
est spectral resolution. The nominal resolving power for 
extended sources is 222, 198, 186 and 175 km/s for the 
[Nelll], H 2 , [Fell] and [SIII] lines, respectively. The ionic 
lines are all resolved and significantly broader than H 2 
which is narrow and unresolved. This result agrees with 
the ground based velocity maps (Fig. ||) and indicate that 
the ionized lines are produced downstream of the shock 
front, while H 2 arises from the precursor (see text, Sect. |J). 




-60 



Fig. 7. Line images of H 2 (2.121 /im) and [Fell] (1.644 
fim) reconstructed from the IRSPEC spectra described 
in Sect. |[ Coordinates are arcsec from the 'ionized 
peak' where all the ISO spectra presented here were cen- 
tered. The sizes of the black rectangles correspond to 
the smallest and largest slits of SWS. The ground based 
[Fell] Al. 644 flux is 70, 82, and 120 (xl0~ 20 W cm" 2 ) in 
the 14x20, 14x27 and 20x33 sqarcsec SWS apertures, re- 
spectively. These numbers are used to compute the aper- 
ture correction factors listed in Table [l| 

of the observed flux of [FeII]A26.0. An alternative expla- 
nation is that the different [Fell] densities simply reflect 
relatively small uncertainties in the atomic parameters of 
Fe + which, for example, overestimate the critical density 
of [Fell] A 1.53 and of the other satellite lines in the near 
infrared. 
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Table 2. Diagnostic line ratios 
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Fig. 8. Energy level diagram showing the paths of the 
[Fell] lines detected by ISO. The near IR lines come from 
higher energy levels which are visualized in Fig. 6b of 
OMD90. 



t 1 ' Log of observed ratio, .including anerture and reddening 
corrections (cf. Table w and Sect. □) 

^ Deduced gas density and temperature, see also Fig. 



From [NeIIi]A3869/H/3=0.41 (Leibowitz & Danziger 1983) 

and [SIII]A9531/H/3=0.16 (Dennefeld |l986j ) adopting the 

case-B H/3 /Bra= 13 ratio 

[SII] lines from Leibowitz & Danziger ( 1983] ) 

Flux of [Fell] 1.644 from IRSPEC image (cf. caption of 

Fig. corrected for extinction using Ah=0.75 



Near infrared [Fell] lines from OMD9C 
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Fig. 9. Contour plots of the [SII] and [Fell] diagnostic line 
ratios listed in Table ||. The curves for [Fell] are com puted 
using the collision strengths of Pradhan & Zhang ( 1993 ) 
and th e transition proba bilitie s of Nussbaumer & Storey 
(|1988|) and Quinet et al. (|l996|) . 



It should be noted that much larger discrepancies 
between the [Fell] densities from FIR and NIR line ratios 
were found from observations of the Galactic center (GC, 
Lutz et al. 1996[ ) where, however, all the [Fell] lines 
but one are reasonably well explained by a collisionally 



excited plasma with T e ~8000 K and 



=10 



3-4 



The only line which does not fit is A 17. 9 which is a 
factor of ~10 lower than predicted. This cannot be easily 
explained neither by alternative excitation mechanisms 
such as e.g. UV pumping (Lucy private communication) 
nor by advocating large uncertainties on the collision 
strengths which are ruled out by the results presented 
here. A partial explanation to the "GC anomaly" could 
be to assume that the foreground extinction at 17.9 /jm 
is a factor of 2-3 larger than that predicted from the 
standard reddening curve. 

The [Nelll], [NeV] and [SIII] ground state lines ra- 
tios are compatible, within the errors, with their low den- 
sity limits, but useful estimates of temperature can be ob- 
tained including measurements of optical transitions from 
the same ions. The values listed in Table || are derived from 
the available fluxes of optical and FIR lines relative to H/3 
and Bra, respectively, and assume a generous error of ±0.2 
dex to include e.g. possible uncertainties in the reddening 
correction of the optical data. The resulting temperatures 
are quite well constrained, nevertheless, and the values are 
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Table 3. Relative abundances 



Lines used' 1 ' 



Abundance ratio' 2 - 1 Solar value 



[FeIII]22.9/[FeII]26.0 


[Fe ++ /Fe + ] = 


= -0.99 






[IN eillj 10.0/ [INeilJ Iz.o 


[In e / in e j 


— — U.4/ 






[ArIII]8.99/[ArII]6.98 


[Ar++/Ar+] = 


= -0.72 






[SIV]10.5/[SIII]18.7 


[S +3 /S++] = 


-1.5: 






[NIII]57.3/[OIII]51.7 


[N++/0++] = 


= -0.6: 


[N/O] = - 


-0.89 


[NeII]12.8/Bra 


[Ne+/H+] = 


-3.7: 


[Ne/H] = 


-3.93 


[Sill] 34.8/[NeII] 12.8 


[Si+/Ne+] = 


-0.66 


[Si/Ne] = 


-0.51 


[PII]32.8/[NeII]12.8 


[P+/Ne+] = 


-2.6: 


[P/Ne] = 


-2.6 


[SIII] 18.7/ [Nelll] 15.6 


[S++/Ne++] 


= -1.0 


[S/Ne] = - 


-0.86 


[ClII]14.4/[NeII]12.8 


[Cl+/Ne+] = 


-2.59 


[Cl/Ne] = 


-2.5 


[ArII]6.98/[NeII]12.8 


[Ar+/Ne+] = 


-1.39 


[Ar/Ne] = 


-1.47 


[FeII]17.9/[NeII]12.8 


[Fe+/Ne+] = 


-0.62 


[Fe/Nc] = 


-0.55 


[NiII]6.64/[NeII]12.8 


[Ni+/Ne+] = 


-1.03 


[Ni/Ne] = 


-1.82 



The most remarkable result is that, apart for the 
anomaly of nickel which is discussed below, all the metal 
abundances are close to the solar values. This indicates 
that the emitting gas is indeed ISM material in which 
the dust grains have been destroyed by the shock front 
thus returning all the refractory species (Si, Fe) into the 
gas phase. It should also be noted that, to the best of 
our knowledge, the ISO data provide for the first time a 
measurement of the abundances of phosphorus and chlo- 
rine in any supernova remnant. Moreover, the abundance 
of silicon from the FIR [Sill] line is much less uncertain 
than those normally derived from UV lines whose intensi- 
ties strongly depend on the assumed gas temperature (e.g. 
Russel & Dopita [l99C| ). 



^ Fluxes from Table |l| including aperture and reddening 
corrections 

Logaritmic values, typical uncertainties are ±0.12 dex 
except for those marked with a column which are ±0.25 dex 



compatible with those predicted for the post-shock region 
of shocks with velocities >300 km/s (cf. Tables 3A-3D of 
DS||). Note in particular that slower shocks are predicted 
to have a post-shock [Nelll] temperature of 15-30 xlO 3 
K and significantly hotter than the observed value. 

The temperature of [NeV] is of particular interest to 
verify if this ion is produced downstream of the shock, 
in which case NeV is predicted to be collisionally ion- 
ized in plasma at about 3.5 10 5 K, or in the much 
cooler (< 2 10 4 K) photoionized precursor. In the first 
case one expects [NeV]A3426/A14.3=8 or, equivalently, 
[NeV]A3426/H/3=0.1, while values a factor of 6-7 lower 
are expected if the emission is from the precursor. Un- 
fortunately, no observation of [NeV]A3426 is available in 
the literature and, given the quite high extinction, mea- 
surements at the required depth (a few % of H/3) are not 
straightforward . 



3.2. Abundances 

A specific advantage of FIR lines is that their emissivities 
depend very little on the gas temperature and abundances 
derived from these lines are, therefore, little affected by 
uncertainties in the assumed value of T e . Also, to mini- 
mize the effect of density variations, we selected for each 
ion the line with the largest critical density which, for 
most species, was 3> 10 3 cm -3 and larger than the values 
of n e derived above. The derived abundances are listed 
in Table [| where the quoted errors also include generous 
uncertainties on temperature (i.e. T e =5000-10000 K and 
T e =7000-15000 K for singly and double ionized species, 
respectively) and density (n e = — 4000 cm -3 ). 



3.2.1. The nickel anomaly 

The relatively large strength of the optical [NiII]A7379 line 
in SNRs and other nebulae has long puzzled astronomers 
because it requires a Ni + /Fe + abundance ratio a factor 
~10 larger than the cosmic value. A similar overabun- 
dance, i.e. [Ni + /Fe + ]=— 0.41 dex or a factor of 7 above 
the solar [Ni/Fe]=— 1.27 value, is inferred here from the 
[NiII]A6.64 ground state transition (cf. Table ||). This re- 
sult is very difficult to understand because the line emit- 
ting gas is shock excited interstellar material which should 
not have, therefore, an anomalous Ni abundance. A num- 
ber of explanations for the optical result have been pro- 



posed (cf. Bautista et al. 1996 for a recent review). These 
include: 

- Selective dust depletion of Fe relative to Ni, which can- 
not hold in our case because most of the iron is in the gas 
phase. 

- Low Fe + /Fe relative abundances, which can be excluded 
because most of the iron is in the form of Fe + (cf. Table ||) . 

- Contribution by UV fluorescence to the [Nill] lines, 
which could be important for the optical lines but has 
negligible effect on the [NiII] A6.64 ground state transi- 
tion (cf. Table 1 of Lucy 1995 ). Also, no obvious source of 
strong U V radiation exists in RCW103 (cf. Sect. 4.1.3 of 
|OMD9C| ). 

- Density stratification with low n e regions emitting [SII] , 
higher density gas emitting [Fell] and the densest (~10 6 
cm~ 3 ) regions dominating [Nill] whose lines have the high- 
est critical density (Bautista et al. |1996| ). However, such a 
high density component should also (but does not) show 
up in other lines with high critical density, in particular 
[NeII]A12.8 which has n cr it = Q 10 5 cm -3 and only a factor 
of 10 lower than the critical density of the [NiII]A6.64 line. 
Since the [Nill] /[Nell] ratio also requires a nickel over- 
abundance of a factor of 6 (cf . Table [3|) , and considering 
also the relatively low densities derived from the [Fell] 
lines (cf. Table ||), one is forced to assume that most of 
the nickel is in regions with densities <; 10 7 cm~ 3 while all 
the other species, including Fe, must be in regions with 
densities <10 cm -3 , which is implausible. 
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We therefore conclude that the apparently large 
Ni + /Fe + relative abundance simply reflects uncertainties 
in the collision strengths for [Nill] whose computed values 
are systematically a factor «10 lower than those of [Fell]. 
In particular, the collision strength of the [NiII]A6.64 
ground state line is only T=0.15 and by far (i.e. a 
factor of >5) lower than the collision strengths of the 
main ground state transitions of any other astrophysical 
abundant species. Adopting, for example, T=1.0 for the 
[NiII]A6.64 transition would yield [Ni+/Ne+]=-1.85 and 
[Ni + /Fe + ]=— 1.18, i.e. values within 0.1 dex of the solar 
Ni/Ne and Ni/Fe relative abundances. 



4. Discussion 

4-1- Comparison with shock models 

Explicit predictions for a few FIR lines have been included 
in the models for relatively slow (<150 km/s) shock inter- 
acting with low density material by Raymond (1979) and 
Shull & McKee ( |l979| ). The first lists [Nell], [Sill] and 
[Fell] while the latter include [Sill], [SIII] and [SIV]. The 
predicted line ratios from the above mentioned models are 
not in good agreement with our results. In particular, the 
observed [SIV] /[SIII] ratio is a factor >5 larger than the 
computed values. Also, most models predict a factor of 
>3 too strong [Sill] (relative to [Nell], [Fell], [SIII]), but 
this could be attributed to uncertainties in the atomic 
parameters of Sill which have been updated several times 
since the publication of the shock model results. The most 
important discrepancy, however, is that the predicted sur- 
face brightnesses are always a factor >10 lower than the 
observed values, which simply reflects the fact that the 
shock of RCW103 is much faster than the values used in 
the above models (cf. the Introduction and below). 

Models of slow shocks interacting with very dense gas 
(n > 10 3 cm~ 3 ) were developed by Hollenbach & McKee 



( |1989| ) who also include explicit predictions for all the 
FIR lines of singly ionized and neutral species. The main 
problem with these models is that they span preshock 
densities much larger than the <300 cm~ 3 required 
to account for the measured electron densities in the 
post-shock region. Consequently, lines with low critical 
densities, e.g. [Sill], are predicted too faint. Moreover, the 
models underestimate the flux of [Fell] lines by a factor 
of -5. 



The most recent models of DS96, which cover shock 
velocities of 200-500 km/s and are more representative 
of the conditions of RCW103, do not however give ex- 
plicit predictions for the FIR lines. Nevertheless, reason- 
ably accurate line ratios can be computed from the values 
of ionic column densities and mean temperatures listed in 
the above paper. To a first approximation, the ratio of two 
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Fig. 10. Behaviour of velocity and precursor sensitive line 
ratios. Th e theo retical values are com pute d from the shock 
models of DS9(: as descri bed in Sect. 4.1. B=Z?/y / n is the 
magnetic parameter (cf. DS96) and the dashed regions 
show the observed values. 



lines from the post-shock region is 
h N 2]2 (T 2 ) 



(1) 



where N and T are the column densities and temperature 
of the emitting ions and j is the line emission coefficient 
which, for FIR lines, is very little dependent on the gas 
temperature. The contribution from the photoionized pre- 
cursor can be also computed from the published tables of 
N and T in the pre-shock region, while the compression 
factor (i.e. the ratio between the electron density in the 
post and pre shock regions) can be estimated imposing 
that the re gions h ave similar H/3 surface brightnesses (cf. 
Sect. 3.2 of |DS96D , i.e. 



\n P N n ] — \n P N n ] . , , 

L P ^precursor L e ylpost — shock 



(2) 



Out of the many FIR line ratios we have identified 
those which are most sensitive to the shock speed and to 
the presence of the photoionized precursor. The behaviour 



of the selected line ratios is plotted in Fig. 10 where the 
most remarkable result is that models including the emis- 
sion from the precursor largely overpredict the strength 
of [SIV] and [NeV]. The velocity dependence of the 
[NeV]/[NeIII] ratio in the post-shock region may appear 
at first sight surprising, but can be easily understood as 
follows. The post-shock [NeV] emission always occurs 
in collisionally ionized gas at T ~ 3.5 10 5 K and whose 
column density, which primarily depends on the shape 
of the gas cooling curve, does not strongly vary with the 



Oliva et al.: ISO infrared spectroscopy of fast shocks. I. Ionized lines in RCW103 



9 



shock speed. The [Nelll] line, on the contrary, could be 
strongly enhanced by emission from photoionized gas in 
the post-shock region, but this only occurs at n > 300 
km/s while slower shocks do not produce enough ionizing 
photons to support a large Nelll zone. In short, the sharp 
decrease of [NeV]/[NeIII] between 200 and 300 km/s is 
because the [Nelll] line emission rapidly increases in this 
velocity interval. 

The main conclusion of this analysis is that all ionic 
lines can be reasonably well reproduced by post-shock 
emission. This conclusion also agre es w ith the broad 
line profiles observed by SWS (Sect. 2.1) and imaging- 
spectroscopy observations of [OIIIJA5007 which show com- 
plex dynamical structures, similar to those seen in [Fell] 
(Fig. |^) and incompatible with emission from the precur- 
sor (Moorwood et al. 1987 ). The only ionic line which 
cannot be accounted for by post-shock emission is [SIV] 
which should be a factor of >10 fainter, but can be re- 
produced adding a quite 'incomplete precursor', i.e. ~5% 
the SI V column density of the precursor predicted by the 
DS96| models (cf. Fig. M). 



Table 4. Comparison between RCW103 (SNR), the 
Galactic center and the Circinus galaxy 



Given the importance that the photoionized precursor 
may have in modelling the spectra of active galaxy nu- 



clei (Sect. 4.3), it is of interest to investigate why little or 
no ionic line emission is observed from the precursor of 
RCW103. We envisage the following possibilities. 

- The precursor in RCW103 is very thin to UV ionizing 
photons, but this is very difficult to reconcile with the fact 
that H2 emission is observed from pre-shock molecular gas 
lying outside of the shock front (cf. Fig. ^). For H2 to ex- 
ist, the molecules must be shielded from the strong field 
of UV ionizing radiation from the shock front or, equiva- 
lently, the preshock region must be optically thick. 

- The shock front in RCW103 is significantly slower than 
so far assumed and below ~150 km/s, the minimum speed 
required to produce a prominent phototoionizcd precur- 
sor. This is in strict constrast with the observed line 
widths and filament dynamics (cf. Sect. § and Fig. |J). 
Moreover, slow shocks cannot explain the very large sur- 
face brightness of the lines which require a large mechan- 
ical power of the shock, i.e. a large nv^ product, n being 
the preshock density and v s the shock speed. More specif- 
ically, the average surface brightness of Bra within the 
ISO beam corresponds to E(H/3)=410 -3 erg cm~ 2 s _1 
sr _1 which, coupled to the predicted values from shock 



models (Eq. 3.4 of pS96| ), yields 

(: 



100 cm- 3 / V 100 km /s 



2.4 



68 cos 6 



(3) 



or, equivalently, a shock speed of about 370 (cos 9) 0A2 
km/s for a pre-shock density of 300 cm -3 . Larger pre- 
shock densities arc effectively excluded by the measured 
electron densities (Table ||) in the post-shock region, i.e. 
after the gas has been compressed by the shock front. The 
factor cos 9 takes into account projection effects such as 



Line ratio 


RCW103 (1) 


GC< 2 > 


Circinus 


[FeIII]A22.9/[FeII]A26.0 


0.089 


3.7 




[Neffl]A15.6/[NeIi]A12.8 


0.89 


0.044 


0.46 


[NeV]A14.3/[Neffl]A15.6 


0.012 




1.0 


[NeIII]A15.6/[FeII]A26.0 


3.1 


1.2 


5.4 


[OiV]A25.9/[FeII]A26.0 


0.30 


0.24 


8.8 


[SiIX]A3.94/[SiII]A34.8 


<0.003 




0.025 



^' Line fluxes from this p aper 
(2) Data from Lutz et al. ( |l996| ) 



(3) Data from Moorwood et al. (|1996; ) 



those modelled in details by Hester (1987) who interpreted 
the bright filaments in IC443 and Cygnus-Loop in therms 
of relatively slow shocks seen quasi edge-on and found 
that small filaments amplified by a factor 10-100 should 
be quite common. However, this model cannot hold for 
RCW103 for the following reasons. This remnant is much 
brighter (factor of >10) than IC443 and Cygnus-Loop. 
The average surface brightness within the relatively large 
ISO-SWS beam (i.e. the value used in Eq. 3) is already 
a factor of ~4 lower than that observed on arcsec scales 
in optical/IR line images of RCW103. The most largely 
amplified edge-on filaments should have smal l radi al ve- 
locities (FWHM<40 km, cf. Fig. 2 of Hester |l987|) amd 
this is not compatible with the observed line widths and 
dynamics. 

- The shock models largely overpredict the contribution 
of the photoionized precursor. Indeed, DS96 state that the 
column density of ionized gas in the precursor might be 
overestimated due to a possibly incorrect treatment of the 
transfer of the UV ionizing photons (cf. end of Sect. 4.2 of 
DS96| ). Moreover, the ionization structure of the precur- 
sor could be much different than computed in DS96 if the 
shock evolves on time scales shorter than ~100 yr, i.e. the 
recombination time in the pre-shock gas. 

4-2. Comparison with the Galactic center 

The region on the line of sight of the GC has a rich spec- 
trum of prominent IR lines which are believed to arise 
from gas with an unusually large Fe gas phase abundance 
and which is primarily photoionized by quite hot stars 
(Lutz et al. 1996 ). Table [| is a comparison between the 
most significant line ratios measured in RCW103 and in 
the GC. 

The [FeIII]A22.9/[FeII]A26.0 ratio is much higher 
(a factor of 42) in the GC spectrum. This implies 
Fe++/Fc+ > 1 and a factor >10 larger than in RCW103, 
regardless of the assumed gas density in the GC. This sim- 
ply reflects the fact that a region predominantly photoion- 
ized by stars, such as those near to the GC, contains only 
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a relatively small fraction of partially ionized gas. The 
recombining region behind the SNR shock front, on the 
contrary, has a large zone of partially ionized gas, which 
is heated by photoionization from the shock front radia- 
tion, and where most of iron is forced into Fe + by the very 
rapid charge exchange reactions with neutral hydrogen. 

The [OIV]A25.9/[FeII]A26.0 ratio is the same in the 
two objects, within the errors. Given the difficulties to 
produce both Fell and OIV with photoionization from 
normal stars, it seems not unreasonable to conclude that 
both species are primarily produced by shock excited gas 
in the line of sight of the GC. 

The [NeIII]A15.6/[NeII]A12.8 ratio is a factor of 20 
lower in the GC than in RCW103 which indicates that 
fast shocks are more effective than late O stars in produc- 
ing Nelll. Moreover, the [NcIII]A15.6/[FcII] A26.0 ratio is 
only a factor of 2.6 higher in the GC than in RCW103 and 
this indicates that a non negligible fraction of the [Nelll] 
emission from the GC could come from shock excited gas. 

4-3. Photoionized precursor and shocks in active galaxy 
nuclei 



According to the shock models of pS9q , the precursor 
could be an important source of lines from high ionization 
species (e.g. [OIII]A5007), but its importance relative to 
the post-shock region may strongly depend on the column 
density of the pre-shock material. In a paper specifically 
dedicated to study the spectral sign atures of shocks in 
active galaxies, Dopita & Sutherland ( 1995| ) consider the 
following limiting cases: 

- Shock only, in which the precursor is very thin and its 
emission is effectively negligible relative to the post-shock 
region. This can fairly well reproduce the line ratios ob- 
served in low excitation AGNs (LINERS). 

- Shock + precursor, where the pre-shock region is opaque 
to the ionizing photons from the shock front. Since the 
ionizing spectrum is quite hard and effectively similar to 
a typical AGN spectrum, the ionization structure of the 
precursor is similar to that of standard narrow line re- 
gions photoionized by the AGN. Consequently, the emerg- 
ing line spectrum is similar to that of standard photoion- 
ization models and could explain, therefore, the high ex- 
citation lines from e.g. type 2 Seyferts. 

In view of this proposed scenario, it is interesting 
to compare the spectra of RCW103 with that of the 
Circinus galaxy, an archetype Seyfert 2 galaxy whose 
observed line ratios are listed in Table [|. The most 
striking difference is that the high excitation (coronal) 
lines are much stronger in Circinus with, in particular, 
[NeV]/[NeIII]=l and roughly 2 orders of magnitude larger 
than in RCW103. Such a strong [NeV] could be in prin- 
ciple compatible with emission from the precursor of a 
u~500 km/s shock (cf. Fig. |l^), while even higher veloc- 
ities, i.e. ^blOOO km/s, could probably account for highest 
ionization coronal lines (e.g. [SiX]). The main problem is 



that these shocks should also emit prominent low exci- 
tation lines from their fast moving post-shock gas, but 
this is incompatible with the observed line profiles which 
are remarkably narrow (FWHM<150 km/s, Oliva et al. 
1994) and similar for all ionization species. Therefore, a 
shock dominated model for the Circinus galaxy seems very 
unlikely and, more generally, the role played by the pho- 
toionized precursor in Seyferts could be questioned on the 
basis of the following arguments. 

If dominated by photoionization, the low excitation 
lines from the post-shock region (e.g. [SII] ) should be 
broader than those from the photoionized precursor (e.g. 
[OIII]), but this is in strict contrast with the observations 
which show that [OIII] and higher excitation lines are 
usually broader than those of [SII] and lower excitation 
species. 

The ISM medium of Seyfert galaxies is well known to 
be quite "porous" , especially within the ionization cones, 
and several arguments indicate that the line emitting 
cloud s are probably density bounded (e.g. Binette et al. 
1996 ). The host galaxies of LINERS, on the contrary, are 
often very rich in both gas and dust, a sp ectacu lar exam- 
ple being NGC4945 (e.g. Moorwood et al. |l996|) . It seems 
therefore curious that the shocks in Seyferts should pri- 
marily impact onto the relatively few large clouds (i.e. 
those with large enough column density to absorb all the 
ionizing photons from the shock) while, in LINERS, the 
shocks should selectively avoid the largest clouds and only 
hit regions with low column densities (i.e. those which can- 
not produce a bright precursor) . 

The absence of significant emission from the pre-shock 
region in RCW103 indicates that shock models may over- 
estimate the importance of the precursor region. 

5. Conclusions 

We have presented IR spectroscopic observations of 
RCW103, a relatively young an d fas t (about 10 3 yr 

supernova rem- 



1984) 



and ~1200 km/s, Nugent et al. 
nant whose prominent line emitting filaments result from 
the interaction of the fast SNR blast wave with ISM 
clouds. The secondary shocks driven into the higher den- 
sity medium have velocities ^300 km/s, as indicated by 
the observed dynamics of the filaments (cf. Fig.||), by the 
very high surface brightness of the lines and by velocity 



sensitive line ratios such as [NeV] / [Nelll] (Sect. 4.1 



The spectrum is dominated by prominent lines from 
low excitation species which have been used to estimate 
the density, temperature and abundances of the emitting 
region. The results indicate this is post-shock gas with 
relatively low density (n e ~10 3 cm -3 ), normal ISM abun- 
dance and essentially free of dust. Although the [Nill] 
lines seem to require a large nickel overabundance we ar- 
gue, however, that this simply reflects uncertainties in the 
atomic parameters and propose that the collision strength 
for [NiII]A6.64 is T=l and a factor of 7 larger than the 



Oliva et al.: ISO infrared spectroscopy of fast shocks. I. Ionized lines in RCW103 



11 



computed value. Smaller uncertainties in some of the 
atomic parameters for [Fell] are also suggested by the data 
which nevertheless exclude large errors in the computed 
collision strengths. 

An important discrepancy with fast shock models is 
that the prominent high excitation lines predicted from 
the photoionized precursor are not seen in RCW103. Since 
low column densities in the pre-shock region are ruled out 
by the observed H 2 emission from the shock precursor, this 
discrepancy probably indicates that shock models largely 
overestimate the strength of the lines from the precursor. 
This may cause problems to shock models of Seyfert galax- 
ies where the high excitation lines should be dominated by 
emission from the photoionized pre-shock region. 
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